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Abstract 

We show that by combining high precision measurements of the atmospheric 5rrF in both the 
electron and muon neutrino (or anti-neutrino) disappearance channels one can determine the neu¬ 
trino mass hierarchy. The required precision is a very challenging fraction of one per cent for both 
measurements. At even higher precision, sensitivity to the cosine of the CP violating phase is also 
possible. This method for determining the mass hierarchy of the neutrino sector does not depend 
on matter effects. 
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Neutrino flavor transitions have been observed in atmospheric, solar, reactor and accel¬ 
erator neutrino experiments. Transitions for at least two different E/L’s (neutrino energy 
divided by baseline) are seen. To explain these transitions, extensions to the Standard 
Model of particle physics are required. The simplest and most widely accepted extension is 
to allow the neutrinos to have masses and mixings, similar to the quark sector, then these 
flavor transitions can be explained by neutrino oscillations. 

This picture of neutrino masses and mixings has recently come into sharper focus with 
the latest salt data presented by the SNO collaboration^. When combined with the latest 
KamLAND experiment P] and other solar neutrino experiments^, 0] the range of allowed 
values for the solar mass squared difference, and the mixing angle, 612 , are ^ 


+ 7.3 X 10“W < 5ml^ < +9.0 x 10“^ 
0.25 < sin^^is < 0.37 


( 1 ) 


at the 90 % confidence level. Maximal mixing, sin^ 0 i 2 = 0.5, has been ruled out at greater 
than 5 a. The solar neutrino data is consistent with ^ and/or z/,-. 

The atmospheric neutrino data from SuperKamiokande has changed only slight in the last 


; 0 | 


are consistent 
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few years j6| and the latest results from the K2K long baseline experiment 
with SK. The range of allowed values for the atmospheric mass squared difference, 6 m ^2 
the mixing angle, 6 * 23 , are 


1.5 X 10 ^eV^ < \Sml 2 \ < 3.4 x 10 +V 


0.36 < sin^ 023 < 0.64 


( 2 ) 


at the 90 % conhdence level. The atmospheric data is consistent with z/^ —>■ oscillations 
and the sign of is unknown. This sign is positive (negative) if the doublet of neutrino 
mass eigenstates, 1 and 2 , which are responsible for the solar neutrino oscillations have a 
smaller (larger) mass than the 3rd mass eigenstate. This is the mass hierarchy question. 

The best constraint on the involvement of the z/g at the atmospheric comes from the 
Chooz reactor experiment and this puts a limit on the mixing angle associated with these 
oscillations, 6 * 13 , reported as 

0 < sin^ 013 < 0.04 (3) 

^ We use the notation of ref.Q throughout. 
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at the 90 % confidence level at = 2.5 x 10“^eV^. This constraint depends on the precise 
value of with a stronger (weaker) constraint at higher (lower) allowed values of 

So far the inclusion of genuine three flavor effects has not been important because these 
effects are controlled by the two small parameters 


6 m? 


21 


Sml^ 


0.03 and/or sin^ 6*13 < 0.04. 


( 4 ) 


However as the accuracy of the neutrino data improves it will become inevitable to take into 
account genuine three flavor effects including CP and T violation. 

One of the goals of the next generation neutrino experiments is to establish the atmo¬ 
spheric mass hierarchy. Many authors have studied how to exploit matter effects in future 
conventional long baseline experiments in supernova explosions Q] or in experiments 
using non conventional neutrino beams produced in a muon collider facility to unravel 
the mass hierarchy. Here we discuss how to make this determination using precision disap¬ 
pearance experiments. 

Genuine three generation effects make the effective atmospheric neutrino measured 
by disappearance experiences, in principle, flavor dependent even in vacuum and thus sen¬ 
sitive to the mass hierarchy and even to the CP phase. This observation suggests an al¬ 
ternative way to access the mass hierarchy by comparing precisely measured values for the 
atmospheric in Ug —»■ (reactor) and —>■ (accelerator) modes. To illuminate this 
rather interesting but experimentally challenging possibility is the purpose of this paper. A 
variant of this idea, using the solar 6 m? scale, can be found in ref. [l^ . 

Assuming three active neutrinos only, the survival probability for the a-fiavor neutrino, 
in vacuum, is given by 


P(n„ ^ n„) = P(z/„ ^ z/„) = 1 - 4|Po3HP„ipsin^ A 

- 4|P„3r|f/a2psin'A 
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( 5 ) 


- 4|P„2r|t/aipsin'A2i, 


where Ajj = 6m‘^jL/4E, 6mfj = mj — m^ and Uai are elements of the MNS mixing matrix, 
jl^ . The three Aij are not independent since the 6 mjj^s satisfy the constraint, 6 ml-^ = 
6m\2 + 6m\i. 

If we define an effective atmospheric mass squared difference, 6 m‘^, which depends linearly 
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on the parameter rj, as follows 


= Sml^ — 7] 5mli = 5ml2 + (1 — h) <5^21 

Sm?L 

so that = A31 - PA21 = A32 + (1 - r])A2i = '1 , (6) 

then we can rewrite Eqn.j^ using the independent variables, A^ and A21, as 

= 4|f/«3p(l - |t/'„3p) [sin^ A^ 

+{ri sin^(pA2i) + r2 sin^((l - r;) A21)} cos 2 A^ 

+ ^{risin(2pA2i) - r2sin(2(l - ?7)A2i)} sin 2 A,, 

+ 4|P„2r|Pairsin2A2i, (7) 


where 


ri 



|U„i| 

2 

|U„i| 

2 +|f /„2 


2 


and r 2 


If/« 2 P 

|77mP + |f/„ 2 P 


1 - ri. 


( 8 ) 


Notice that the coefficient in front of sin2A^ is the derivative of the coefficient in front of 
cos2A^, with respect to r]A 2 i, up to a constant factor. Therefore by choosing rj so as to 
set the coefficient in front of sin 2A,, to zero one also minimizes the coefficient in front of 
cos2A^. That is, if rj satishes 


T] = 


2A 


arctan 


21 


r 2 sin 2 A 


21 


A, 


(9) 


ri + r 2 cos 2 A 21 

one minimizes the effects of both sin 2A,^ and cos 2A^ terms and this Sm^ with rj k, r 2 is 
truly the effective atmospheric 5rn?, 6mlfi\a, measured in Ua disappearance experiments. 
The approximation rj = r 2 is excellent provided that A 21 1. 

Using this approximate solution for t], the effective atmospheric for the a-flavor is^ 




|t/m| 

P + |[/„2| 

^ 5m^2 

|Um| 

P + Ua2 

2 


= ri 5m\^ + r2 


( 11 ) 


^ An alternative way to derive this is to notice that the first extremum, of the terms in Eqn. proportional 
to |C/a 3 p, occurs when 


lUiipAai + |[/Q,2pA32 _ TT 

|C/alP + |C/„ 2 P " 2 ’ 


( 10 ) 


to first non-trivial order in A 21 . 
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then the full neutrino survival probability in vacuum, Eqn[S], can be rewritten as 

1 - P{Ua Va) = 4|17„3p(l - It/asH ^eff 

+{ri sin^(r2A2i) + r2 sin^(riA2i)} cos2Aefr 
-^{r 2 sin( 2 riA 2 i) - n sin( 2 r 2 A 2 i)} sin 2Aeff 
+ 4|f/Q2nf^aiP sin^ A 21 . ( 12 ) 

If the coefficients in front of the cos 2Aeff and sin 2Aefr terms are expanded in powers of A 21 , 
one hnds 

{risin^(r 2 A 2 i)+r 2 sin^(riA 2 i)} = rir 2 A 2 i + C>(A 2 i) 

^{r2 sin(2riA2i) - ri sin(2r2A2i)} = ^rir2(r2 - ri)A^i + 0(A^J, (13) 

and one can see clearly that all terms linear in A 21 have been absorbed into the Agg terms. 
This conhrms that Eqn^J, is the effective atmospheric 5m^ to hrst non-trivial order 

in Note also that the hrst term odd in Agg occurs with a coefficient proportional to 

A 21 which, at the hrst extremum, is a suppression factor of order 10“'^. 

To understand the physical meaning of the ehective atmospheric 5m? it is useful to write 
it as follows 


5m^g|„ = ml - (m^)i 2 . 


where 


\U^2\ 

+ \Ual\ 

pmi 

\u^i\ 

p + 

\U^2\ 

2 


(14) 


Now (m^)i 2 has a clear interpretation, it is the a-havor weighted average mass square of 
neutrino states 1 and 2. Thus the ehective atmospheric 5m? is the diherence in the mass 
squared of the state 3 and this havor average mass square of states 1 and 2 and is clearly 
havor dependent. 

The three havor average mass squares are^ 


(me)i2 = ^["^2 + ^ 1 - cos 26 'i 2 5 ml^] 

= -[^2 + + (cos26*12 — 2COS5sin013sin26*12tan6*23) 5^21] (15) 

{ml) 12 = - [ml + m\ + (cos 20i2 + 2 cos 5 sin 0i3 sin 20i2 cot ^ 23 ) 5 ml]\, 


^ Dropping terms of order sin^ 9iz8rn?i 
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± 5m^ ± 5m^ ± 5m^ 



FIG. 1: The vacuum survival probability, P{i'a as a function of EjL for the two mass 

hierarchies using three different choices of the atmospheric whose flips sign, with constant 
magnitude, changes the hierarchy: 5rnl^\a (left panel), (middle panel) and 5m\2 (right panel). 
The survival probability for the two different hierarchies coincide to high precision when the effec¬ 
tive 5m? ’s, Ean|TCl ITT] . are used (left panel) whereas they differ appreciably with the other two 
definitions. For this figure we have used sin^ 023 = 0.5 (maximal mixing), sin^ 0 i 3 = 0.04 (Chooz 
bound), sin^ 012 = 0.31, 5m‘^i = -1-8.0 x 10“^ eV^ and the atmospheric 5m? to be 2.5 x 10“^ eV^. 

where the r-flavor flavor average is given for completeness only. 

It is now obvious that and disappearance experiments measure different In 

fact the three are^ 

Smlgle = cos^ 012(5^3]^-|-sin^ 0i25m32 (16) 

= sin^0i2(5m33-|-cos^ 0125^32-|-cos(5sin013 sin 2012 tan0235’^2 i 

Smlff I T = sin^ di 2 Sml^ + cos^ 0 i 25 m 32 — cos 5 sin 0 i 3 sin 20 i 2 cot 023(5m2i. (18) 

^ The effective atmospheric mass squared difference for the muon channel has been discussed in ref. [idf. 
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FIG. 2: The fractional difference of the electron and muon neutrino effective atmospheric <5m^, 
Ae^ = (|(5mgg|e — |(5mgg|^)/|(5mgg|, as a function of sin^ ^13 for the normal and inverted hierarchies 
showing the dependence on cos <5. The vertical scale varies linearly with the not so well known 
ratio of 5m\i/5rn^2^ here we have used 6 m\i = 8.0 x 10“^ eV^ and 5 m^2 — 2-5 x 10“^ eV^. In a 
reactor i/g disappearance experiment, precision measurement of the effective atmospheric 
is probably very difficult unless sin^ @13 > 0.005. 

In Fig. [T] we show the survival probability in the z/g and disappearance channels using 
three different choices of the atmospheric 5m?‘ whose sign flip, with constant magnitude, 
changes the hierarchy from normal to inverted. When we use for the a flavor, the 

change in the survival probability is very small when we flip the hierarchy i.e. the magnitude 
of this Smljj is insensitive to which hierarchy nature has chosen. Although (<^’^ 32 ) 

works better for z/g (z/^) disappearance experiments neither choice is as good as 
Thus, in summary, Equ^Hl; is the atmospheric measured by Ue disappearance 

experiments and Eqn[T7j, is the atmospheric measured by z/^ disappearance 

experiments upto corrections of 0 { 6 mli/ 6 ml 2 )^- 

Whether the absolute value of is larger or smaller than the absolute value of 
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depends on whether is larger or smaller than The relative magnitude of these 

two Srn^ is determined by whether the mass squared of the 3-state is larger or smaller than 
the mass squared of the 1- and 2-states, i.e. by the neutrino mass hierarchy. It is easy to 
show that the difference in the absolute value of the e-flavor and /i-flavor is given by 

= ±(5m2j^(cos 20i2 — cos (5 sin 0i3 sin 20i2 tan 023); (19) 


where the -|- sign (— sign) is for the normal (inverted) hierarchy. Thus by precision mea¬ 
surements of both of these one can determine the hierarchy and possibly even cos 6 at 
very high precision. This identity, Eqn. in. is the principal observation of this paper. 

In Fig. 121 we show the fractional difference in the effective atmospheric for the normal 
and inverted hierarchy, as a function of sin^0i3. For the normal hierarchy, independently 
of 6, this normalized ratio is always positive, while for the inverted hierarchy, it is always 
negative. While the size of difference between the two hierarchies is smallest for cos 5 = 1, 
for this value of S, the difference between the two hierarchies increases as sin^ 0i3 goes to 
zero, as can be seen from FaufTH]. 

What kind of precision is required? Given that 


6mn 


21 


5ml2 


— and 
30 


cos 2012 ~ 0.38, 


( 20 ) 


the difference in the magnitude of the two effective atmospheric is 1 to 2%. Currently, 
the uncertainty on the size of this difference is dominated by the experimental uncertainty 
on the ratio of the solar to atmospheric 5m?'’s. To determine the hierarchy we need to 
determine whether |5mgg|e is larger, normal hierarchy, or smaller, inverted hierarchy, than 
I^TTiggl^. Thus determining the hierarchy with a conhdence level near 90% one needs to 
measure both to better than one per cent precision. These are very challenging levels 
of precision for atmospheric measurements both within a given experiment and between 
two different experiments. In Fig|21 we have calculated the required precision as function 
of the C.L., measured in sigmas, assuming that the two experiments have the same % 
precision. From this hgure we see that for a 90% C.L. determination of the hierarchy one 
would require ~0.5% precision on both 5mlj-f measurements. Achieving such precision will 
require signihcant innovation. 

So far our discussion has only been in vacuum. What about matter effects? How much 
do they shift the hrst extrema? For the z/g disappearance channel the shift in the extrema 
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FIG. 3: The required percentage precision need to determine the neutrino mass hierarchy verses 
the confidence level of that determination. Here we have assumed both effective atmospheric 6m? 
are measured with the same precision, Uee = The cosine of the CP violating phase is varied 
from +1 (bottom) through 0 (dashed line) to -1 (top). Again, the vertical scale varies linearly 
with the not so well known ratio of 5m2i/5m‘^2- For this figure we have used 0.032, the same as in 
Fig.m 

is proportional to (aL) where a = CpNe/'^^ ~ (4000 km)“^. Thus the expected shift is less 
than 0.1% for a baseline of a few kilometers. The size of this shift has been conhrmed by 
a numerical calculation. For the disappearance channel again the shift in the extrema 
is again proportional to (aL) but here the baseline could go up to 1000 km. However the 
coefficient in front of (aL) is proportional to sin^ 2^13 and cos 26 * 23 /cos^ ^23 both of which 
are small numbers. Using an energy so that the hrst minimum occurs at 1000 km, we have 
calculate numerically the size of the shift assuming sin^ 26*13 is at the Chooz bound and 
found that the maximum shift is 0.4%. This maximum shift occurs when 6*23 is as larger as 
is allowed by atmospheric neutrino data. If sin^ 6*23 and/or sin^ 6*13 are smaller than these 
maximum values then the shift is smaller. Also the shift at baselines smaller than 1000 km 


9 











are proportionally smaller. Therefore, we conclude that in general matter effects can be 
safely ignored, or corrected for, in disappearance experiments whose baseline is less than 
1000 km. 


In summary we have demonstrated that high precision measurements of the effective 
atmospheric in both the 9e (reactor) and —>■ (long baseline accelerator) 

channels can determine the neutrino mass hierarchy independent of matter effects. The sign 
of the difference determines the hierarchy. For any reasonable conhdence level determination 
the precision required in both channels is a very chal enging fraction of 1%. The next 
generation of long baseline experiments such as T2K Q and NOz/A 0 estimate their 
precision on the effective atmospheric Sm^ at 2%. However, so far there has been no physics 
reason to push this to a precision measurement. For the reactor channel the emphasis so 
far has been on the observation of non-zero 613 very little effort has been made on 
a precision determination of the effective atmospheric This kind of precision, can 

perhaps be achieved in beta beam facility Q| . We realize that to make these measurements 
to the precision suggested is very challenging experimentally. However we encourage our 
experimental colleagues to give this some thought especially since this method has a different 
dependence on the unknown CP violating phase, cos 5 versus sin 5, compared with long 
baseline experiments. 

While we were completing this manuscript, ref.^^ appeared which discusses the physics 
of this possibility in a pure 3-flavor frame work as well as discussing other possible ways of 
determining the hierarchy. 
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